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ABSTRACT 

o 

Linear eigenmodes of a spherically symmetrical ultra-relativistic blast wave (the 
| Blandford-McKee, BMK, solution) are calculated. The BMK solution is shown to be 

q ' stable and strongly non-universal. It is stable because all the eigenmodes decay. Non- 

universality of the BMK solution (BMK is not an intermediate asymptotic) is a conse- 
quence of causality. In terms of eigenmodes - some eigenmodes decay too slowly. For 
each degree I, there exists an eigenmode which decays at the smallest rate. The am- 
plitude of this eigenmode, defined as relative perturbation of energy behind the shock 
front, is constant at early time. Later, when the blast wave slows down to T ~ 1/10, an 
oscillatory decay commences, and the amplitude drops to less than 10% of the initial 
value at T ~ Z/100. The non-universality is surprisingly strong. Near the end of the 
ultra-relativistic stage, perturbations with I <^ 200, that is more than 10,000 differ- 
ent harmonics, are suppressed by less than a factor of 10. Spherical symmetry is only 
reached by the time the blast wave slows down to non-relativistic velocities, when the 
^ c*j Sedov/Taylor/von Neumann solution sets in. 
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Subject headings: shock waves - gamma rays: bursts 



1. Ultra-relativistic blast waves, GRBs, and their aftreglows 

It is thought that ultra-relativistic blast waves occur naturally, and GRB afterglows come from 
synchrotron-emitting ultra-relativistic blast waves (Piran 1999). 

A spherically symmetrical ultra-relativistic blast wave (the BMK solution) was given by Bland- 
ford and McKee (1976). We will show that this solution is stable but strongly non-universal (§3). 
In §2 we derive the basic equations, it includes a simple derivation of the BMK solution. 

Qualitatively, non-universality follows from causality. For the BMK solution, the shock front 
moves with Lorentz factor T oc t~ 3 / 2 . Consider a light signal propagating along the front, starting 
at a polar angle 8 = 0. Then, with c = 1, (dr/dt) 2 + (rd6/dt) 2 = 1, which, for r > 1, gives 
8 = |r _1 . Two regions of the shock front, separated by angle 6, do not talk to each other until the 
blast wave slows down to T ~ O^ 1 . 

It turns out that non-universality is very strong (§3). The non-universality of ultra-relativistic 
blast waves should be observed as non-universality of early (ultra-relativistic) GRB aftreglows. 
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2. Basic equations. The BMK solution. 

Relativistic hydrodynamics equations are written as d k T ik = (Landau & Lifshitz 1987) . 
Here x k = (t,r) are Cartesian coordinates tied to the rest frame of the unshocked fluid. The 
energy-momentum is T lk = (Au l u k — g lk )p in the shocked fluid. The energy-momentum in the 
unshocked fluid, T tk , has only one non-zero component, and without loss of generality in the final 
result, we take r 00 = 1. The shock front is at F = 0, the boundary conditions follow from the 
hydrodynamics equations: 

T ik d k F = r ik d k F, F = 0. (1) 

The 4- velocity is u k = 7(1, v). In spherical polar coordinates, without loss of generality in the 
final result, v = vf + uQ. We use the three-dimensional form of the hydrodynamics equations: 

d T 00 + 8 a T 0a = 0, (2) 

and 

d T 0a + dpT* = 0, (3) 
which we project onto f and 9. Only first-order terms in u, dgp, dgv are kept. We get 

3t[(4 7 2 - l)p] + d r [4~f 2 vp] + 8r _1 7 2 «P + 4 7 2 p(r sin0) -1 flk[sin0u] = 0, (4) 

d t [^ 2 vp\ +d r [(4j 2 v 2 + l)p] + 8r~ 1 j 2 v 2 p + 4~f 2 vp(r smO^delsmOu] = 0, (5) 
dtl^pu] + d r [4j 2 vpu] + 12r~ lr y 2 vpu + r~ 1 dop = 0. (6) 

We parameterize the shock front as F = r — R(t, 9) = 0. Then eq. (1) gives, up to first order 
in u, dgR: 

A 1 2 p{d t R-v]-pdtR = d t R (7) 

A-i 2 v[d t R-v] - 1 = 0. (8) 

A-f 2 u[d t R -v] + R^deR = 0. (9) 

We introduce a new independent variable, x = t — r (dt — > dt + d x , d r — > —d x ). We will assume 
x <C i, r. We parameterize the shock position as R = t — x s (t, 6). Denote 7 2 = q. Then, 

» = 1 -^-8? + 0( «" 3) (10> 
Keeping two leading orders in g -1 , we get 

dt[(4q- \)p] +d x [(l + —)p] +8r- 1 (q- -)p + Aqp(r sine)' 1 d e [smeu} = 0, (11) 
at[(4g - 2)p] + - i-)p] + Sr- 1 ^ - l)p + 4(g - ^)p(r sinflr^sinflu] = 0, (12) 



d t [4qpu] + d x [(2 + i-)p«] + ^r" 1 ^ - ^)p« + = 0. 

In the leading order in g -1 , these can be written as 

dt(qp) + 2t~ 1 qp+ -^d x p = -qp(t sin0) -1 d0[sin0u], 

dtp + At^p + dxiq^p) = -2p{ts\n6y 1 d e {sm6u} 1 

1 u 1 

d t u + r x u + — d x u + - — d x p + - — -d e p = 0. 
2g 4gp 4qpt 



The boundary conditions, at x = x s : 



1 1 d e x s 



Ad t x s 3d t x s 



2.1. BMK solution 



Under spherical symmetry, u = 0, <9# = 0, eqs.(14), (15), and (17) are written as 

<%(<?p) + -qp+\d x p = 0, 

%> + jP + d x {q~ x p) = 0, 



Ad t x s 3d t x s 
These have a self-similar BMK solution 

_ £2 _ 2r 2 _ 

9 ~2? P " 3^712' Xs "8f2' 



x s ' 8nt 3 ' 



E is the energy of the explosion. 



3. Linear eigenmodes 

We write linear perturbations of the BMK solution as 

^ = q(x,t)Y l0 {6), S -P=p( x , t )Y l0 (9), 5 -^ = S(t)Y m (9), 5u^u(x,t)^^. 
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We use x, and t as independent variables (<% — > dt — 4t 1 x^%' ^ — ► 8r 2 i 1 9 x ). From (14)-(16), 
we get for the perturbations 

td t q + td t p-4 X d x q+¥- q = l(l+l)u, (24) 

tfltP + 12 X 5 x p - 16 X 9 x g + f q = 21(1 + l)u, (25) 
td t u + 4x9 x u - H u = _ (26) 
The boundary conditions at % = 1 are 

q = -jtd t 5, p = -jtd t S+^5, u = gf2- ( 27 ) 

We perform the last change of independent variables, t — > e*, % — > e 4x . We denote i and x 
derivatives by an over-dot and a prime. With no loss of generality in the final answer, we take such 
an explosion energy that T = t~ 3 / 2 . The final form of equations and boundary conditions is: 

x>0: q + 3q' - 3p' -q = -1(1 + l)u, (28) 
p + Sp' - Aq' + fq = 21(1 + l)u, (29) 
i + u '_^ u = _i e 3*+4xp. (30) 

x = 0: p = q-\q, u= ^e 3t (p- q). (31) 

There are two different cases. Spherical, I = 0, and non-spherical, I ^ 0. All non-spherical 
perturbations are isomorphic, because (28)-(31) are invariant under /(/ + 1) — > Xl(l + 1), u — > u/A, 
*-»t-lnA/3. 



3.1. Spherical perturbations 

For spherical perturbations, eqs. (28)-(31) have a time independent solution, q = 0, p = 1, 
u = 0. This is the slowest decaying eigenmode. This eigenmode is trivial - the difference between 
BMK solutions of different energies. There are two other / = eigenmodes. One can check that 
(28)-(31) have solutions oc e ~ At ~ Ax and oc e ~ 5t ~ 6x . These decay algebraically (that is algebraically 
in the old time variable). These eigenmodes are also trivial - the differences between the BMK 
solution and BMK solutions propagating from t ^ and r^0. Numerical simulations (§3.2) show 
that BMK solution is the intermediate asymptotic of a spherical explosion. 



3.2. Numerical simulations 

Numerical simulations of eqs. (28)-(31) were performed after introducing "Riemann invariants" 
(Landau & Lifshitz 1987). Let 

f = q + fp, 9 = q-^P- (32) 

Then 

/ = (2V3 - 3)/' + (1 - + g ) + (V3 - l)Z(i + IK (33) 

3 = _(2>/3 + %' + (| + j=)(f + g)-(V3 + 1)1(1 + l)u, (34) 

il = - u > + ^ u -J- e ^ U _ g) . (35) 

Since g and u propagate to larger x, one needs two boundary conditions at x = 0; these are 
given by eq. (31). We need one boundary condition at the end of the simulation domain x — x max . 
We checked that this boundary condition is irrelevant if x max is large enough. 

Results, the time evolution of p(x = 0,t), are shown in fig. 1. 
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Fig. 1. — Relative pressure perturbation (normalized to relative pressure perturbation at t = 0) as 
a function of 1(1 + 1) /T 2 . 



